Notebook criado para obter a temperatura central antes e depois dos "pulsos termicos". Esta informação será utilizada em um parametro do MESA que modifica a quantidade de camadas em relação a temperatura central
In [1]:
# lendo os dados do mesa
import mesa as ms
import numpy as np
folder2 = '/home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2/'
s = ms.history_data(folder2)
In [2]:
s.cols
Out[2]:
In [4]:
# extraindo os modelos, temperatura central e quantidade de zonas
model = s.get('model_number')
zones = s.get('num_zones')
ct = s.get('log_center_T')
In [24]:
# grafico de zonas vs model (para confirmar os dados de outro notebook)
import matplotlib.pyplot as plt
%matplotlib inline
plt.rcParams["figure.figsize"] = (15., 10.)
#fig = plt.figure(figsize=(15., 10.))
plt.plot(model, zones, 'b.')
plt.plot(model[499], zones[499], 'ro')
plt.plot(model[2089], zones[2089], 'go')
print 'camadas modelo 500: ', zones[499]
print 'camadas modelo 2090: ', zones[2089]
O grafico acima esta de acordo com os resultados obtidos no notebook "quantidade de camadas"
In [9]:
# grafico temperatura central vs model
fig = plt.figure(figsize=(15., 10.))
plt.plot(model, ct, 'b-')
plt.plot(model[499], ct[499], 'ro')
plt.plot(model[2089], ct[2089], 'go')
plt.ylabel(r'$\log T_c$', fontsize = 16)
plt.xlabel('Model', fontsize = 16)
print 'log temperatura central modelo 500: ', ct[499]
print 'log temperatura central modelo 2090: ', ct[2089]
O Kepler pediu para obter a temperatura central no meio do caminho das supergigantes para modificar a quantidade de camadas a partir desse ponto. Determinando o modelo pelo HRD:
In [22]:
# plot do HRD
t = s.get('log_Teff')
L = s.get('log_L')
fig = plt.figure(figsize=(15., 10.))
plt.plot(t, L)
plt.plot(t[154], L[154], 'r.')
plt.plot(t[2089], L[2089], 'go')
plt.gca().invert_xaxis()
plt.title(r'HRD $\alpha_{MLT}=2$', fontsize=18)
plt.xlabel(r'$\log T_{Eff}$', fontsize=16)
plt.ylabel(r'$\log L$', fontsize=16)
Out[22]:
In [28]:
# determinando a temperatura central do modelo 155
plt.plot(model, ct, 'b-')
plt.plot(model[154], ct[154], 'ro')
plt.plot(model[2089], ct[2089], 'go')
plt.plot([0, 3000], [8.6, 8.6], 'k--')
plt.ylabel(r'$\log T_c$', fontsize = 16)
plt.xlabel('Model', fontsize = 16)
print 'log temperatura central modelo 155: ', ct[154]
print 'log temperatura central modelo 2090: ', ct[2089]
Pelos resultados acima, a faixa de temperatura dos "pulsos termicos" eh de 8 a 8.6 (log T_c).
Com isso, vou definir os parametros:
Pelo o que eu entendo, nessa faixa de temperatura o parametro mesh_delta_coeff sera dividido por tres, triplicando o numero de camadas e com sorte obtendo os pulsos termicos. Vou rodar o teste no pc da UFRGS e salvar na pasta alpha2_3x_camadas
Ideia: Ver se tem como manter registro do numero de camadas enquanto roda o ./rn . Talves os parametros num_trace_history_values = 1 e trace_history_value_name(1) = 'num_zones' de certo
In [ ]: